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Abstract 



> 

O , 

i An asymptotic stability analysis of spatially homogeneous models of 

f*"^ . Bianchi type containing tilted perfect fluids is performed. Using the 

known attractors for the non-tilted Bianchi type universes, we check whether 
' they are stable against perturbations with respect to tilted perfect fluids. 

We perform the analysis for all Bianchi class B models and the Bianchi 
type VIo model. In particular, we find that none of the non-tilted equilib- 
, rium points are stable against tilted perfect fluids stiffer than radiation. 

O" 1 We also indicate parts of the phase space where new tilted exact solutions 

^' ( ' might be found. 



1 Introduction 



'_ ' In order to explain the isotropy and homogeneity of the universe and investi- 



gate the form of its small expansion anisotropics, it is important to understand 
the complete spectrum of anisotropic cosmological models. The subset of the 
most general anisotropic and inhomogeneous cosmological models that are most 
amenable to a full analysis are the spatially homogeneous Bianchi universes 
because the Einstein field equations for these cosmologies reduce to ordinary 
differential equations. They have been studied in considerable detail by many 
authors but attention has largely focussed on the behaviour of models with 
comoving perfect fluids. In this paper we will extend these analyses to the situ- 
ation where the fluid is 'tilted' with respect to the normals to the hypersurfaces 
of homogeneity and the fluid is not comoving. 

In the study of Bianchi type universes there is special interest in the par- 
ticular exact solutions of the Einstein equations which are attractors for more 
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general classes of solution (see P for a review). The two recent papers, |21E], 
have completed the investigation of all non-chaotic Bianchi models with a non- 
tilted 7-law perfect fluid and their asymptotic attractors are understood. This 
understanding takes the form of a phase-space description of the dynamics of 
the autonomous ordinary differential equations describing the cosmological evo- 
lution. However, when we turn to Bianchi models with tilted fluid motion 0], 
the situation is far from completely understood. Up to three additional degrees 
of freedom can be added to the phase spaces in these cases. Some analyses have 
been made of Bianchi universes with tilted fluids pUH| , in particular type 77 [7| 
and type V [HIEIEIj but a more general analysis has not been by date. The 
object of this paper is to provide a local stability analysis of a wide range of 
Bianchi universes containing tilted perfect fluids in order to elucidate their late- 
time behaviours. Specifically, we will consider Bianchi class B models, together 
with the type VIq model of class A, filled with a tilted 7-law perfect fluid. We 
will base our analysis on the results from the non-tilted analysis. The main 
reason for this is that, due to the complexity of the equations of motion, a com- 
plete analysis would need to find, or rule out the existence, of new tilted exact 
solutions which might act as attractors. However, a relatively simple analysis 
based on the results from the non-tilted models allows us to go a significant way 
towards understanding the behaviour of tilted universes. 

Tilted universes display many features that are of considerable interest in 
the study of Einstein's equations and their physical implication. They include 
rotating universes and are therefore of importance for any attempts to formulate 
and evaluate Mach's Principle ^5> the conditions needed for the appearance of 
closed timelike curves, and the properties of quantum cosmologies that tunnel 
from 'nothing'. The changes in the cosmological evolution of anisotropies that 
arise in brane- worlds are also of special interest but, so far, the influence of 
non-comoving velocities have not been studied. They are likely to be the most 
important factor in the evolution as can be seen from the effects of anisotropic 
pressures studied by Barrow and Maartens [T2). 

Before we embark on the formal analysis of the full Einstein equations gov- 
erning the tilted Bianchi universes it is helpful to use some elementary physical 
arguments to determine when we might expect critical changes in the behaviour 
of velocities and vorticities to occur as the equation of state varies. Using a 
fairly simple argument we can give a rough estimate of the stability properties 
to be expected of non-comoving perfect fluids A perfect fluid with equation 
of state, p = (7 — l)p, and constant 7, will evolve according to p oc a~ 37 cx 9 2 , 
where a(t) is the mean scale factor and 9 = 3a/ a is the volume expansion rate. 
Consider an expanding universe that is close to isotropy: an eddy of fluid, with 
mass m, will conserve angular momentum, Ilu, if 

Ma 2 uj = constant. 

For tilted fluids with velocity v — auj and mass m oc pa 3 , this implies 

pa 5 ui = pa A v = constant. 
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Thus for perfect fluids 



v oc a 37 " 4 , u oc a 37 " 5 . 



The ratio of the distorting rotational energy density uj 2 to the isotropic density 
evolves as 

oc — oco a7 iu . 
Z p 

We see that we expect particular values of 7 to mark thresholds in the evolution 
of models containing non-comoving velocities and vorticity. In particular a 
threshold exists when 7 = 4/3: rotational velocities grow with a(t) for 7 > 
4/3 |13|. The influence of vorticity on the expansion dynamics grows when 
7 > 10/9. The significance of the value 7 = 10/9 can be seen, for example, 
in the exceptional model, discussed section 13.51 In the analyses to follow we 
find that these critical values of 7 recur, along with the value 7 = 2/3, which 
marks a transition between universes that isotropise at late times (for 7 < 2/3) 
and those that do not (for 7 > 2/3). Of course, the detailed behaviour is more 
complicated than this discussion implies because the Bianchi types need not 
be close to isotropy and there are geometrical constraints on how the different 
components of the 3- velocity are able to evolve if energy and momentum are to 
be conserved . 

In the next section we will determine the equations of motion for the tilted 
fluid in Bianchi universes. Then, in scction|21 we will perform a stability analysis 
for the various asymptotic behaviours that are known from the study of non- 
tilted universes. In section ^ we will analyse the type Vlh model. Finally, in 
section we summarise and discuss our results. 



2 The equations of motion 

We choose a time variable so that the time is measured along time- like geodesies 
orthogonal to the surfaces of transitivity. The surfaces of transitivity are there- 
fore spatial and we say that the spacetime is spatially homogeneous. For the 
Bianchi cosmologies - which admit a simply transitive symmetry group acting 
on the spatial hypersurfaces - we can always write the line-element as 

ds 2 = -dt 2 +g ab {t)u) a u) h , 

where g ab is a non-singular symmetric 3x3 matrix only dependent upon time, 
oj a is a triad of one-forms obeying 

1 This simple physical argument can be extended to the early stages of brane- world cos- 
mologies. In brane-worlds we have 6 oc p rather than 9 2 oc p at early times and so the 
distortion to the metric created by any vortical motions is much more significant than in gen- 
eral relativistic cosmologies as the universe expands. We have w/H oc a 67-5 and rotational 
distortions grow in all brane-worlds with 7 > 5/6 and this includes the cases of matter and 
radiation domination. 
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and C a bc are the structure constants of the Bianchi group type under consid- 
eration. The structure constants C a bc can be split into a vector part a b , and a 
trace-free part n ab by 

C a bc = e bcd n da - 5\a c + 6 a c a b . 

The matrix n ab is symmetric, and, using the Jacobi identity, a b = (l/2)C a ba is 
in the kernel of n ab 

n ab a b = 0. 

The Bianchi types are divided into two classes: 

• Class A: a c = 0. 

• Class B: a c 7^ 0. 

First, we will concentrate on the class B spacetimes. In this case we can 
choose the orientation of the triad so that a c = aS 1 c . This implies that n ab is of 
rank 2 and can be written 

" 

n 22 n 23 . 

n 23 n 33 

Consider a perfect fluid with energy-momentum tensor 

= (p + p)Ufj,u v + pg^. 

The normalised 4-velocity of the fluid m m is in general tilted with respect to 
the normals to the hypersurfaces of homogeneity, n a — t a ; and the fluid three- 
velocity v a , where w M = (u°, v a ), will be assumed to be non-zero. The equations 
of motion for the fluid, T^ u = 0, can be written as a closed system of equations 

-{p+p)V a V a = pv a V a (1) 

-(p + p)(u Q v c + v a v b C a bc ) = pu°v c (2) 
(p + p)(u° + v b a b + u°9) +u°p = 0. (3) 

Here, 9 is the volume expansion factor 9 — g/2g. In our analysis, eqs. (lll.'il 
are the essential equations because we will consider the stability of non-tilted 
solutions, that obey the Einstein equations with v a = 0, with respect to the 
perturbations introduced by the present of non-comoving velocities. 

Recently, the two remaining Bianchi models with non-tilted 7-law perfect 
fluid have been analysed HH2|; so an non-tilted Bianchi models have now been 
analysed and their asymptotic behaviours determined. Using these results, 
which give all the attractors in the set of non-tilted perfect fluid cosmologies, 
we can investigate the stability of the late-time non-tilted attractors against 
perturbations from tilted matter. 
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It is useful to introduce a orthonormal triad e a a (t) by defining 



9ab, 




We can now write eqs. Ill HI) in terms of the components «a = e a a v a in the 
orthonormal frame: 



The stability of a particular solution can now be checked by linearising these 
equations around the exact non-tilted solution describing the asymptote in ques- 
tion. For non-tilted equilibrium points, the linearised version of eqs. il llol) will be 
a closed set of equations, and hence, the stability analysis can to a large extent 
be done by considering these equations alone (provided the stability analysis 
has already been performed for the remaining equations). 

More precisely, our stability analysis is based on the following observation. 
The total dynamical system (including the Einstein equations) can be written 
close to an equilibrium point, say at x = 0, 



where F(x) is at least quadratic in x. In our case, when considering a non-tilted 
equilibrium point, the matrix A will be of the block form 



This means that the eigenvalues of A are the union of the eigenvalues of Ai and 
A2. Thus if we choose equilibrium points that are stable in the set of non-tilted 
cosmologies, we only have to check the eigenvalues for Ai; i.e. the linearised 
version of cq. (j5j- This simplifies the stability analysis considerably and turns 
out to be very fruitful. However, there are some limits to this procedure: 

1 . If the real parts of one or more eigenvalues of Ai are zero, then the stability 
analysis is inconclusive from analysing Ai alone. 

2. For tilted equilibrium points, A will in general not have the form J7J); i.e. 
the complete set of equations is needed to establish its stability. 

Keeping these limitations in mind we shall perform a perturbation analysis 
for various spacetimes. 




x = Ax + F(x), 
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3 Stability analysis 



We will study the various class B cases, plus the type VIq case which is closely 
related. We will henceforth assume that the tilted perfect fluid obeys the 
barotropic equation of state 2 

V = (7 - 

Note that ever-expanding cosmologies with inflationary types of fluids with < 
7 < 2/3, will all asymptote to a FRW model with a(t) cx t 2 / 3j as t — > 00 
in accord with the cosmological "no-hair" theorem ^5]^. These fluids will 
therefore not be considered explicitly in what follows. However, the 7 = 2/3 
will appear as a stability boundary for the anisotropic attractors as expected. 

The notion of tilt does not apply to an exact vacuum solution. Hence, for 
equilibrium points corresponding to vacuum spacetimes we can only define a tilt 
when they are perturbed away from any exact vacuum solutions. Thus, in what 
follows, when refer to stability with respect to the introduction of tilt we are 
determining whether the universe at late times is best described by a tilted or a 
non-tilted model. Also, when referring to the vacuum solutions we will always 
mean the non-tilted equilibrium points. 



3.1 Bianchi type VII h 

The most general Bianchi universes containing the open Friedmann universe as a 
subcase are those of type VIIh ■ The general late-time asymptotes for non-tilted 
spacetimes of type VIIh are the plane- wave vacuum solutions of Lukash [201 - 3 
For these solutions [23 



= dx 



u> = e sx [coswrrdy — sinwxdz] 
a; 3 = e sx [sin toxdy + cos uixdz] 
s(l-s) = w 2 sinh 2 2/3, 



(8) 
(9) 
(10) 
(11) 



and 



e" = 



t 

t s e f3 cos (w \nt) -i s e^sin(a;lnt) 
t s e~ 13 sm{u hit) t s e - ' 9 cos(wlnt) 



The group parameter, h, is given by 



2 Non-perfect fluid stresses, for example anisotropic pressures, can have significant effects 
on the evolution of anisotropics; see for example refs. II 711 21 . 

3 These vacuum solutions can easily be generalised to plane- wave solutions with an electro- 
mag netic field 
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Furthermore, 9 = (2s + l)/t, and they reduce to isotropy when s = 1 and 
L) = j3 = 0. Eq. (JSJ) is to lowest order in v& 

p + j8p = 0. (12) 

This is the same as for non-tilted fluid. Defining expansion-normalised energy- 
density by SI = p/0 2 , we require ft < for the solutions to be future stable. 
This yields the bound 

2 

This is the usual bound required for the stability of these solutions in the space of 
non-tilted Vllh universes with perfect fluid found earlier [Tl ll8l l^]. Linearising 
eq. (J5J, using eq. I|12|l . and dividing by jp, we get the following system of 
equations for the 3- velocities: 

v[ = [(2s + 1)7 - (2s + 2)] vi 

^ = [(2a + 1)7 - (3« + 1)] «j - e^owg 

4 = [(2s + 1) 7 - (3s + 1)] v s + e" 2 ^. (13) 

Here, we have also introduced a new time- variable t — hit and a prime denotes 
derivative with respect to r. The eigenvalues of this system of equations, which 
determine the stability of the plane-wave solutions, are 

Ai = (2s + 1) 7 - (2s + 2), 

A 2 , 3 = (2s + l)7-(3s+l)±tw. (14) 

Hence, the plane-wave solutions of type Vllh are stable when 

2 3s + 1 



< 7 < 



2s + 1 2s + 1 

Note that (3s + l)/(2s + 1) < 4/3; so the solutions are unstable for fluids stiffer 
than radiation. This reflects the tendency for very stiff fluids to 'spin up' as 
they expand in order to conserve angular momentum |13j. as discussed above. 

For radiation, 7 = 4/3, the solutions are unstable for s < 1. The case 
s = 1 yields a zero eigenvalue for radiation and has to be considered separately. 
However, s = 1 implies that the solutions are LRS and are therefore the same as 
the LRS type V solutions. Type V evolution has been considered elsewhere, and, 
in particular, the stability of this solution has been analysed in detail |51l31ll(J| . 
For generic perturbations 4 this solution is found to be unstable. Hence, all 
plane- wave solutions of type Vllh are unstable for radiation (7 = 4/3). 

Since for a given s, any value of h is possible, we can conclude that for any 
given h there are plane-wave solutions of type Vllh that are stable in the set 
of Vllh models, whenever 

2 4 
- < 7 < -. 

3 ' 3 

A summary of the results for type Vllh is illustrated in Fig. ^ 

4 In I1UI . it is shown that this solution is, in fact, an attractor for a set of non-zero measure 
of solution trajectories in the state space of irrotational type V models. However, for some 
perturbations it is unstable and is therefore generically unstable. 
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Figure 1: Stability diagram for the Bianchi type VII^ Lukash plane- wave so- 
lutions. Here, the labelled curves are defined by the following functions: A is 
7 = 23TT> S is 7 = and C is 7 = |f±f . The region marked "Stable" is 

bounded by the curves A and B, while the region marked "Unstable" is bounded 
by A, B and C as shown. The latter region has two directions in which the so- 
lutions are unstable against tilt. All the remaining unmarked regions of the 
diagram are also those in which the Lukash solutions arc unstable. 



3.2 Bianchi type V 

The general late-time attractor for non-tilted type V models is the isotropic 
Milne universe. This vacuum solution can be extracted from the Lukash plane 
waves above by choosing ui — 0, s — 1. Hence, from the above we see that this 
solution is stable provided that 

2 4 
- < 7 < -. 

3 ' 3 

This agrees with earlier results of other authors 6 8 MTU] . 

3.3 Bianchi type VI h 

Consider the type VIh plane-wave solutions given by 



8(1-8) 



= dx 

= e^+^dy 

= eS s -^ x dz 

= b 2 + r 2 , 



(15) 
(16) 
(17) 
(18) 



and 



e „ = 



t 
t s+b 






_ r j.s— b 
b b 

t s-b 



The group parameter is given by 



and 8 = (2s + l)/t. Following the same procedure as in the type Vllh case, we 
obtain a lower bound from eq. (JSJ): 



l + 2s 



< 7- 



Linearising eq. ||SJ we get 



V 3 



[(2a + l)7-(2a + 2)]w 1> 
[(2s + 1)7 - (3s + 1 + b)] v 2 , 
[(2a + l)7-(3a + l-6)]« a , 



and the eigenvalues are: 



(19) 



(20) 



Ai = (2s + 1) 7 - (2s + 2), 
A 2 , 3 = (2s + 1)7- (3s + l±6). 

Assuming that bound (|19|l holds, we have 



(21) 



A 3 > 2 - (3s+ 1 - b) > -(1 - V^h). 
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Hence, all models with \h\ < 1 are unstable. The threshold model, of LRS type 
III = VI-i, has one solution with a zero eigenvalue. 

Combining these results, we see that for \h\ > 1 there are stable plane wave 
solutions whenever 



2(1 -ft) 1- 
< 7 < 



V^ft- 4ft 4 

1 -3ft < 3' 



1 - 3ft 

Consider next the Collins Vlh perfect fluid solution 
u) 1 = dx 



2 f(2-y)+e 

dy 



e 2 ~< 

r(2--,)-c 



or = e 2 -* x dz 
c 2 = (2- 7 )(3 7 -2), 



and 



where < r < 1, and 



2 — f + i- 

t~^r 



2—1- 



(2- 7 ) 2 r 2 



(22) 

(23) 

(24) 
(25) 



This solution is stable to the future for non-tilted perfect-fluid cosmologies when- 
ever 



2 2(1 - ft) 

- < 7 < — -. 

3 1 1 - 3ft 



Linearising eq. (0, we get 



7-2 



(26) 



We note that 



— (57 - 6 - rc) t>2 

— (57 — 6 + rc) Ua. 
27 6 



Ai < 0, A 2 < A 3 
for 7 < 2 and so stability is assured if 

(5 7 -6) 2 



or in terms of 7, 



A3 < 



7 < 



'^ < (3 7 -2 ) 2 - 

2(3 + ^/^h) 
5+3V^ft ' 



(27) 
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When — h < 1, this bound in tighter than the upper bound in cq. <|26[) : when 
—ft > 1 the bound Il26[l is the tightest. Hence, we can conclude that there are 
future stable Collins type Vlh perfect fluid solutions in the presence of tilt so 
long as 



— < 7 < mm 
3 ' 



2(1 -ft) 2(3 + V-h) 



6 

<5- 



3.4 Bianchi type IV 

The late-time attractors for non-tilted type IV models are also plane waves. 
Their stability can be obtained from the type Vlh analysis by taking the limit 
b — > to obtain the eigenvalues 



Ai = (2s + 1) 7 - (2s + 2), 
A 2 , 3 = (2s + 1) 7 - (3s + 1). 

Hence, there will be stable solutions whenever 



(28) 



- < 7 < -. 

3 1 3 



3.5 The exceptional model of Bianchi type VP 



1/9 



The asymptotic dynamics of the exceptional model has been studied with a 
non-tilted perfect fluid by Hewitt 22 , and more recently by Horwood et al 
A stability analysis was done for both non-tilted fluid and with fluids with a 
single non-zero tilted component by Barrow and Sonoda [H]- 
Consider the Robinson- Trautman vacuum solution 



dx 



2^6 



dz 



(29) 
(30) 
(31) 



and 



2 L 









This solution has unusual volume expansion, with 6 = 9/(5t). This solution is 
an attractor for all non-tilted perfect fluids with 10/9 < 7. This requirement 
for stability is easily seen to come from eq. (j!2|) . Linearising eq. JHJ), we obtain 



14)«j 

io)« 2 

4)« 8 . 



(32) 
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However, due to the exact vanishing of the 02-equation of Einstein's field equa- 
tions, we have to set — 0. Hence, the type VI*_ X , 9 implies that only and 
can be non-zero. Thus we have only two eigenvalues to consider 

Ai = i(9 7 -14), A 3 = |(3 7 -4). 

The first of these eigenvalues was found by Barrow and Sonoda [Hj. However, A3 
gives an even tighter bound for stability. Thus, the Robinson- Trautman solution 
is stable if 

10 4 

9 " ' 3 

For 7 = 10/9, the attractor for the non-tilted case is the Wainwright Vl*_ 
7 = 10/9 perfect fluid solution j2H] which interpolates between the Collins 
perfect-fluid solutions, studied earlier, and the Robinson- Trautman solution. 
The solution is 



and 



u) 1 = dx 

u> = e s x dy 

u> = e 5 dz 





" t 


0" 




wt t~s 







. 


t\ 


2 = 


= 9r 2 /4- 


1. 



eigenvalues 
- 4 A - 2 



(33) 
(34) 
(35) 



Here, we have 9 = 9/ (5t), and w 
Linearising eq. 0), we get the 



Since both are negative, this solution set is asymptotically stable as well. 

When 7 < 10/9, the non-tilted solutions asymptote to the Collins yi_i/g 
perfect-fluid solution studied earlier. However, here we must again exclude 
one of the eigenvalues because of the vanishing of the (02) Einstein equation. 
Nonetheless, all eigenvalues are negative when h = —1/9 and 7 < 10/9, so this 
does not tighten the bounds already arising in the non-tilted analysis. Hence, 
there are stable solutions of the exceptional model when 

2 10 
- < 7 < — . 

3 ' 9 

3.6 Bianchi type VIq 

This model can be obtained as the h — > limit of type Vlh universes, both 
geometrically and dynamically. In the non-tilted case, the Collins VIq perfect- 
fluid solution is the future attractor for 2/3 < 7 < 2. Allowing for tilted fluids 
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we see that the Collins VIq perfect-fluid solution is stable whenever 

2 6 

3 ' 5 



4 Tilted Bianchi type Vlh equilibrium points 

An interesting point emerges from the analysis of the type Vlh models. We can 
see that the instability of the Collins perfect-fluid solution in the exceptional 
model with 7 — 10/9 corresponds to an instability of the Collins solutions with 
respect to a tilted fluid for more general h (see Fig. where this is illustrated by 
the fact that the line B intersects h = —1/9 at 7 = 10/9). Hence, at h = —1/9, 
the instability with respect to tilt is exchanged for an instability of the shear: 
when h = — 1/9 there does not need to be a matter source to drive the solution 
away from the Collins solutions. This may indicate that the exceptional model, 
in the presence of tilted fluid, is not as dynamically exceptional as one might 
think and that the rest of the type Vlh models will behave in a similar way. 
Thus there should exist exact tilted solutions of type Vlh, analogous to the 
Robinson- Trautman vacuum solutions, which act as attractors. Such solutions 
are not known in general, but we expect them to exist. 

Hence, one is tempted to search for solutions on the form 



uj = dx 

a, 2 = e^^dy 
a; 3 = e^+^dz 



(36) 
(37) 
(38) 



and 



e „ = 



* 
w-it t a 
w 2 t t T 



where u>i,2, c, r, s and b are constants. Due to the power-law dependence of 
the triad, this metric is self-similar by construction and the Einstein equations 
reduce to purely algebraic equations. However, the equations are still quite 
complex and difficult to solve in full generality. Rosquist |23|, and Rosquist and 
Jantzen have found some solutions of Bianchi type VIq of this type. In 
general, such spacctimes will describe solutions with tilt and non-zero vorticity. 
Furthermore, due to the complexity of the solutions no analysis of their stability 
has been done. Preliminary results indicate that the solutions found by Rosquist 
and Jantzen may play an important role in the late-time behaviour for at least 
a class of tilted type VIq models |2fi| . Notwithstanding the lack of a complete 
analysis, we believe that these kinds of solutions may be important for the 
late-time behaviour of models in the remaining region in Fig. [21 namely that 
bounded by the curves 7 = 4/3, A, and B. One piece of supporting evidence 
for this conjecture is, as already mentioned, the fact that the exceptional case 
intersects this region. We believe that the Robinson- Trautman solution shares 
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Figure 2: Stability diagram showing the stable solutions for the Bianchi type 
Vlh cases. The labelled curves are defined by the following functions: A is 7 = 



-h-ih 



, Bis 7 



and C is 7 



2(1-M 
l-3/i 1 



The region marked "Vacuum 



l-3/i ' " ' 5+3\/=ft 

plane waves" is bounded by the curves A and C, while the region marked "Collins 
perfect fluid" is bounded by 7 = 2/3, B and C. The region bounded by 7 = 4/3, 
A and £? contains none of the solutions we have investigated as attractors, 
except for the exceptional case which cuts through this region at h = — 1/9. For 
non-tilted cosmologies the region corresponding to "Collins perfect fluid" is the 
entire region bounded by 7 = 2/3 and C; while for "Vacuum plane waves" it is 
the region above C. 
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many common features with these "missing" solutions (apart from the fact that 
this is the only vacuum solution). 

5 Discussion 

A summary of our results are given in table ^ 

We have studied the stability of all the exact non-tilted class B and type VIq 
Bianchi-type universes which act as attractors for the general non-tilted cosmo- 
logical evolution at late times. This has led to a determination of their local 
stability in the presence of non-comoving motions fluid motions. None of the 
non-tilted attractors are stable with respect to perfect fluids stiffer that radia- 
tion (p > p/3). The type Vlh models show a fairly complex dependence on the 
equation of state and the group parameter defining the 3-curvature anisotropy 
of the universe. However, there may be indications from this analysis that there 
are some "missing" stable models of this type which have not yet been found. 
The exceptional FI*]/ g model is the only case with known exact solutions of 
this type. Further analysis of this "missing" set of solutions is required. For 
type Vllh, the situation is simpler: whenever 2/3 < 7 < 4/3 there are always 
some Lukash plane- wave metrics that are stable for any given h =/= 0. 

This investigation has been a step towards acquiring a complete understand- 
ing of all tilted Bianchi models. Clearly, there are some missing pieces in the 
analysis: class A universes require a separate analysis and amongst those of 
class B the most prominent lacuna is the "missing" set of solutions of type Vlh- 
Also, we have only done a local stability analysis of non-tilted equilibrium point. 
The type V analysis JU| shows that the evolution in the whole state space may 
be far more complex than local perturbative analyses can reveal. New equilib- 
rium points with non-zero tilt could exist which act as additional attractors. 
Furthermore, due to the nature of our method, the stability of the solutions 
along the boundary between stable and unstable regions (i.e. those that have 
one or more eigenvalues with zero real parts) remains unsettled. These models 
require further analysis. The Bianchi type III universes appear interesting in 
this context. The type III dust model (7 = 1) is the centre model in Fig. [21 
where the curves A, B and C all intersect. Hence, the neighbouring models 
appear to have a very sensitive dependence on the defining parameters. 

This study is has been a first attempt to extend our systematic understanding 
of the dynamics of Bianchi type universes to those with tilted fluid motions. It 
should be viewed as a first step towards gaining that understanding. We hope 
that the lacunae we have defined and the stability analyses we have presented 
will provide a basis for the development of a complete picture of the behaviour 
of tilted Bianchi type universes. Although non-comoving fluid motions have 
generally been neglected in studies of anisotropic models of the early universe, 
they are by no means negligible I28j . Under a wide range of realistic equations 
of state they have been shown to be the most important source of anisotropy in 
expanding universes. In the case of brane-world cosmologies j^Z] the effects of 
these motions will be even more dominant and no reliable understanding of the 
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Solution Type 


Matter 


Stability 


Vllft (Lukash) 


I < 7 < f 




Stable 
Unstable 


IV (p-wave) 


\ <7<| 
3 < 7 < 2 


Stable 
Unstable 


V (Milne) 


i < 7 < 4 

3 ' 3 

3 S 7 S * 


Stable 

I n CtT" O r~\ 1 

U llbLcLUlc 


Vl h (Collins) 


3 7 <. mm ^ 1 _ 3h , 5+3v/ ^ 7I y 

2(3+v / ^7I) ^ 2(1-/!) 
5+3v^7I ~~ ' ~~ l-3ft 


Stable 
Unstable 


Vlh (p-wave) 


2(l-h) ^ ^ l-A/=B-4ft 
l-3ft " ' " 1-3/t 

i ^lr M <7<2 


Stable 
Unstable 


VP 1/9 (Collins) 


!<7<f 


Stable 


vr 1/9 (w) 


7- ^ 

' 9 


Stable 


VI* 1/9 (R-t) 


9 ^ ' ^ 3 

^<7<2 


Stable 
Unstable 


VI (Collins) 


1 <7< ! 
!<7<2 


Stable 
Unstable 



Table 1: A summary of our results: the exact solutions listed in the first column 
are all attractors in the non-tilted case for the equation of state and model 
parameters indicated in second column. The third column records whether they 
are stable or unstable when tilted fluids are introduced. Here, W=Wainwright, 
and R-T=Robinson-Trautman. 

early evolution of these cosmologies will be possible in general without a full 
analysis of tilted fluid motions and vorticity. 

Acknowledgments 

SH was funded by the Research Council of Norway and an Isaac Newton Stu- 
dentship. 

References 

[1] C.G. Hewitt and J. Wainwright in Dynamical Systems in Cosmology, eds: 
J. Wainwright and G.F.R. Ellis, Cambridge University Press (1997) 

[2] J.T. Horwood, M.J. Hancock, D. The, J. Wainwright, |gr-qc /02 10031 

[3] C.G. Hewitt, J.T. Horwood, J. Wainwright, |gr-qc/021107T| 

[4] A.R. King and G.F.R. Ellis, Gommun. Math. Phys. 31 (1973) 209 



16 



[5] O.I. Bogoyavlenskii and S.P. Novikov, Sel. Math. Sov. 2 (1982) 159; origi- 
nally published as Trudy Sem. Petrovsk. 1 (1975) 7 

[6] J.D. Barrow and D.H. Sonoda, Phys. Reports 139 (1986) 1 

[7] C.G. Hewitt, R. Bridson, J. Wainwright, Gen.Rel.Grav. 33 (2001) 65 

[8] I.S. Shikin, Sov. Phys. JETP 41 (1976) 794 

[9] C.B. Collins, Comm. Math. Phys. 39 (1974) 131 

[10] C.G. Hewitt and J. Wainwright, Phys. Rev. D46 (1992) 4242 

[11] D. Raine and M. Heller, The Science of Space-Time, Pachart Publ. Tucson, 
(1981). 

[12] J.D. Barrow and R. Maartens, Phys. Rev. D59 (1998) 043502, 1 

[13] J.D. Barrow, Mon. Not. R. astron. Soc. 179 (1977) 47P 

[14] G.F.R. Ellis and M.A.H. MacCallum, Comm. Math. Phys. 12 (1969) 108 

[15] M. Dcmianski and L.P. Grishchuk, Comm. Math. Phys. 25 (1972) 233 

[16] M. Demianski, Relativistic Astrophysics, Pergamon Press, Oxford, (1985) 

[17] J.D. Barrow, Phys. Rev. D 55 (1997) 7451 

[18] J.D. Barrow and D.H. Sonoda, Gen. Ret Grav. 17 (1985) 409 

[19] R.M. Wald, Phys. Rev. D28 (1983) 2118 

[20] V. Lukash, Sov. Phys. JETP 40 (1975) 792 

[21] S. Hcrvik, gr-qc/0210080 

[22] C.G. Hewitt, Gen. Rel. Grav 23 (1991) 691 

[23] J. Wainwright, Gen. Rel. Grav. 16 (1984) 657 

[24] K. Rosquist, Phys. Lett. 97A (1983) 145 

[25] K. Rosquist and R.T. Jantzcn, Phys. Lett. 107A (1985) 29 

[26] S. Hervik, A Class of tilted type VI o universes, in preparation. 

[27] J.D. Barrow and R. Maartens, Phys. Lett. 532B (2002) 153 

[28] C. Uggla, H. van Elst, J. Wainwright, and G.F.R. Ellis, gr-qc/0304002 



17 



